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INTRODUCTION 


This  research  is  part  of  the  recent  AFOSR  initiative  in  Biomimetics,  which  has  the  broad 
overall  goal  of  seeking  to  elucidate  and  understand  the  salient  characteristics  of  naturally 
occurring  biological  composite  materials,  with  an  aim  toward  mimicking  these  features  in  devising 
improved  man-made  composite  materials.  The  biological  composite  of  interest  in  this  effort  is 
bone  tissue,  or  more  specifically  Haversian  compact  bone  tissue.  A  major  underlying  task  in  this 
project  is  to  develop  a  more  deterministic  understanding  of  the  relationships  between 
microstructural  details  and  macroscopic  mechanical  behavior.  The  unique  fiber/matrix  interphase 
material,  which  is  distinctive  of  this  type  of  bone  tissue,  is  of  particular  interest  because  of  the 
recognized  significance  of  interfaces  in  developing  new  higher  performance  composites.  The 
approach  combines  composite  micromechanics  modeling  techmques  to  analytically  model 
structure/property  relationships  with  companion  experimental  testing  of  tissue  samples  to  evaluate 
and  validate  the  modeling. 

The  specific  objectives  of  the  project  as  stated  in  the  original  proposal  were: 

♦  To  develop  a  two-dimensional  micromechanics  finite  element  model  for  Haversian 
bone  tissue  that  includes  the  main  microstructural  components  (interstitial  bone, 
cement  lines,  and  secondary  osteons  with  separate  lamellar  bone  layers  and  a 
Haversian  canal); 

♦  To  validate  the  model  by  comparing  predictions  with  established  trends  reported  in 
the  literature  and  by  conducting  a  series  of  experiments  measuring  mechanical 
properties  as  a  function  of  microstructural  characteristics; 

♦  To  use  the  model  to  study,  identify,  and  characterize  basic  structure/property 
relations  exhibited  by  bone  tissue  and  formulate  ways  in  which  this  knowledge  can  be 
applied  to  the  development  of  more  advanced  aerospace  composite  materials;  and 

♦  To  identify  and  specify  further  improvements  and  modifications  needed  in  the  model 
and  additional  experimental  work  required  for  more  definitive  validation  of  the  model. 
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RESEARCH  DESCRIPTION 


Background 

Bone  is  a  composite  material  consisting  of  both  solid  and  fluid  phases.  Blood  and 
extracellular  fluids  comprise  the  fluid  phase.  The  solid  phase  is  composed  primarily  of  a  collagen 
matrix  interspersed  with  a  mineral  compound,  hydroxyapatite.  Bone  tissue  is  generally  divided 
into  two  types:  trabecular,  or  cancellous,  bone  and  compact,  or  cortical,  bone.  Compact  bone  is 
much  more  dense  than  trabecular  bone,  which  possesses  a  spongy  structure.  Three  types  of 
compact  bone  are  generally  recognized:  woven,  lamellar,  and  Haversian  (or  osteonal).  Further 
details  of  these  distinctions  are  well  documented  [1,  2,  3].  Haversian  compact  bone  is  the 
predominant  bone  type  in  the  mid-shaft  region  of  long  bones  like  the  femur  and  tibia  in  adult 
humans.  The  primary  microstructural  element  of  Haversian  bone  is  the  secondary  osteon,  as 
shown  in  Fig.  1.  The  very  first  osteons  formed  in  growing  primary  lamellar  bone  are  termed 
primary  osteons.  Osteons  formed  from  subsequent  remodeling  are  termed  secondary  osteons. 


Figure  1.  Idealized  view  of  bone  mid-shaft  showing  compact  bone  microstructure 
(adapted  from  [4]). 

Osteons  are  basically  cylindrical  in  shape  with  diameters  ranging  from  100|xm  to  200p.m. 
They  are  packed  rather  tightly  together  and  oriented  roughly  parallel  to  one  another  along  the  axis 
of  the  long  bone.  The  sectional  view  of  Fig.  1  exaggerates  the  relative  size  of  osteons  since  there 
are  typically  10  to  50  layers  of  osteons  through  the  wall  thickness  of  the  long  bone.  Osteons  are 
typically  composed  of  20  to  30  concentric  layers  of  lamellar  bone  with  a  hollow  center,  the 
Haversian  canal.  Haversian  canals  provide  passageways  for  nerves  and  arteries.  Volkmann  canals 
run  transverse  to  the  main  bone  axis  periodically  connecting  adjacent  Haversian  canals.  A 
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relatively  thin  (l-5|xm)  layer  of  highly  mineralized  material,  called  the  cement  line,  surrounds  each 
secondary  osteon.  The  cement  line  is  laid  down  as  bone  resorption  ends  and  new  bone  formation 
begins.  This  amorphous  'interphase'  material  also  contains  polysaccharides  and  collagen.  Small 
cavities  called  lacunae,  contain  bone  cells  and  lie  between  adjacent  lamellar  bone  layers.  Lacunae 
are  interconnected  with  one  another  and  to  Haversian  canals  through  tiny  passageways  called 
canaliculae.  The  remaining  space  between  Haversian  systems  is  occupied  by  less  organized 
'interstitial'  lamellar  bone.  In  mature  human  bone,  this  bone  tissue  is  actually  composed  of 
fragments  of  osteons  that  have  been  resorbed  away  (incompletely).  Woven  bone  and  pnmary 
lamellar  bone  can  also  persist  as  interstitial  bone.  Mature  equine  (horse)  bone  is  also 
predominantly  Haversian,  but  bovine  (cow  or  steer)  bone  is  typically  a  mix  of  Haversian  and 
lamellar  bone  (also  called  plexiform  bone). 


The  particular  type  of  bone  tissue  studied  in  this  research  is  Haversian  compact  bone 
tissue  because  its  microstructure  resembles  that  of  a  fiber  reinforced  composite  material.  The 
fibers  in  this  case,  however,  are  Haversian  systems,  or  secondary  osteons.  Thus,  they  are  hollow, 
layered,  and  surrounded  by  a  distinct  'interphase'  fiber/matrix  material,  known  as  the  cement  line. 
The  hollow  fibers  are  furthermore  not  perfectly  straight  but  in  fact  spiral  slightly  and  branch.  The 
'matrix'  material  is  the  interstitial  bone  lying  between  complete  osteons.  The  niatrix  is  therefore 
quite  similar  in  composition  and  properties  to  the  fiber  component.  A  microradiograph  of  equine 
Haversian  compact  bone  tissue  is  shown  in  Fig.  2, 


Figure  2.  Microradiograph  of  Haversian  compact  bone  transverse  section. 
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I 

Theoretical  Modeling 

I  The  underlying  strategy  in  seeking  to  develop  theoretical  nucromechanics  models  relating 

macroscopic  mechanical  behavior  to  microstructural  details  has  been  to  build  upon  techniques 
I  developed  for  aerospace  composites.  The  two  basic  methodologies  examined  are.  (1)  finite 

’  element  unit  cell  analyses,  and  (2)  analytical  micromechanics  approaches  (particularly  effective 

medium  methods). 


Porosity  effects.  A  two-dimensional  finite  element  model  was  developed  based  upon  an 
idealized  representation  of  Haversian  compact  bone  microstructure  [5],  The  bone  microstructure 
is  depicted  in  Fig.  3.  The  three  solid  components  are  the  osteons,  interstitial  bone,  and  cement 
lines,  which  represent  the  fiber,  matrix,  and  interphase  constituents,  respectively.  Note  as  well 
that  the  fiber  is  hollow  with  the  Haversian  canal  being  the  only  void  space  represented.  Fiber 
packing  geometries  have  been  studied  for  either  a  square  array  or  hexagonal  array,  the  latter  of 
which  is  shown  in  Fig.  3.  A  typical  mesh  generated  for  this  case  is  shown  in  Fig.  4.  The  finite 
element  modeling  pre-processor  PATRAN  was  used  in  conjunction  with  the  ABAQUS  finite 


element  solver  to  analyze  this  model. 


Baseline  mechanical  properties  used  for  each  constituent  are  summarized  in  Table  1.  The 
interstitial  bone  is  taken  to  be  slightly  higher  in  modulus  than  a  typical  osteon  because  it  is  older 
than  the  more  recently  formed  osteons  and  therefore  more  highly  mineralized.  Table  2  compares 
the  hexagonal  array  results  with  those  generated  using  a  square  array.  Note  that  En  is  the 
longitudinal,  or  axial,  modulus  whereas  E22  and  E33  are  moduli  in  the  transverse  directions.  The 
Haversian  porosity  in  this  case  is  3%.  Note  that  the  results  are  essentially  the  same  for  hexagonal 
and  square  packing  with  the  exception  of  a  slight  discrepancy  for  E22..  In  addition,  hexagonal 
packing  offers  the  advantage  of  being  able  to  model  higher  fiber  volume  percentages,  so  it  has 
been  chosen  as  the  method  for  all  subsequent  finite  element  unit  cell  modeling. 


Table  1  -  Constituent  properties  used  in  micromechanics  model 


Table  2  -  Predicted  macroscopic  properties  for  hexagonal  vs.  square  packing  geometry 


Note:  elastic  moduli  are  in  GPa 


The  Haversian  porosity  is  a  parameter  which  greatly  influences  the  macroscopic  material 
properties  of  bone  and  is  considered  an  important  variable  in  micromechanics  modeling.  In  the 
finite  element  model,  the  Haversian  porosity  is  altered  by  changing  the  diameter  of  the  Haversian 
canal,  i.e.  the  inner  radius  of  the  hollow  fiber.  Porosity  variations  from  3%  to  9%,  which  are 
reasonable  values  found  in  the  literature,  have  been  studied.  Table  3  details  the  influence  of 
porosity  on  transverse  and  longitudinal  moduli  as  well  as  Poisson's  ratio.  The  results  are 
baasically  as  expected:  The  moduli  generally  decrease  with  increasing  porosity,  and  the  transverse 
moduli  are  more  sensitive  to  porosity  variations. 
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Table  3  -  Effect  of  porosity  on  macroscopic  properties 


Porosity  (%) 

Ell 

E22 

V12 

V23 

V21 

3 

12.3 

11.7 

0.30 

0.30 

0.29 

5 

12.2 

11.1 

0.30 

0.29 

0.28 

7 

11.9 

10.6 

0.30 

0.29 

0.27 

9 

11.6 

10.1 

0.30 

0.28 

0.26 

Note;  elastic  moduli  are  in  GPa 


The  Mori-Tanaka  effective  medium  method  of  estimating  the  overall  mechanical 
properties  of  a  fiber  reinforced  composite  has  been  extended  to  allow  for  porosity  by  Zhao,  et  al. 
[6],  Even  though  experimental  data  is  limited,  they  report  reasonable  predictions  for  variations  in 
mechanical  properties  as  a  function  of  porosity.  The  simplest  case  to  consider  as  a  model  for 
Haversian  compact  bone  is  that  of  'needle-shaped'  voids  representing  the  Haversian  canal.  That 
is,  the  voids  are  long  and  slender,  parallel  to  one  another,  and  aligned  with  the  overall  longitudinal 
aks.  The  variation  in  longitudinal  modulus  with  porosity  is  linear  in  this  case,  however,  whereas 
experimental  data  suggests  a  power  law  relation  [7].  More  plausible  comparisons  with 
experimental  data  are  nevertheless  possible  by  using  the  Mori-Tanaka  method  for  3-D  ramdomly 
oriented  spheroidal  voids  with  appropriate  aspect  ratios.  Results  are  shown  in  Fig.  5  for  sever^ 
aspect  ratios.  The  ratio  of  the  longitudinal  modulus  to  the  modulus  of  the  solid  phase  (En/Eo)  is 
plotted  versus  porosity.  The  aspect  ratio  (a)  of  10  represents  a  prolate  spheroidal  void  and  the 
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other  three  aspect  ratios  (.10,  .05,  .02)  represent  various  oblate  spheroidal  voids.  Prolate  voids 
are  more  'rod-like',  while  the  oblate  voids  are  more  'disk-like'.  The  power-law  expression  used  to 
fit  the  experimental  data  is  given  by: 


E  =  33.88(1-P)‘“®^  (inGPa) 


The  curves  in  Fig.  5  that  best  fit  the  data  are  for  oblate  voids.  This  certainly  does  not  describe  the 
shape  of  the  Haversian  canal,  but  the  experimental  data  is  for  bovine  bone,  which  contains  both 
Haversian  and  plexiform  types  of  compact  bone.  No  experimental  data  was  previously  available 
for  modulus  versus  porosity  for  exclusively  Haversian  compact  bone,  so  the  expenments 
conducted  as  part  of  this  research  have  specifically  targeted  this  need. 

Most  attempts  to  develop  analytical  expressions  relating  porosity  and  mechanical  properties  of 
various  engineering  materials  have  resulted  in  empirical  correlations  based  upon  extensive  experimental 
data.  These  expressions  typically  take  the  form  of  power  law  relations  with  particular  coefficients  and 
exponents  for  different  materials.  In  a  similar  vein,  Bert  [8]  developed  a  semi-empirical  approach 
combining  theoretical  and  empirical  considerations.  Cylindrical  or  spheroidal  voids  can  be  modeled 
with  square  or  hexagonal  array  spacings.  Elastic  modulus  predictions  are  compared  in  Fig.  6  for  the 
finite  element  unit  cell  analysis,  Mori-Tanaka  methods,  and  the  semi-empirical  method  of  Bert.  The 
Mori-Tanaka  predictions  are  for  the  case  of  needle-shaped  voids.  The  moduli  values  are  normalized  by 
dividing  by  the  elastic  modulus  of  the  solid  phase.  All  three  approaches  yield  the  same  results  for  the 
elastic  modulus  in  the  axial  direction,  i.e.  along  the  main  longitudinal  axis  of  the  long  bone.  The 
cylindrical  voids  are  also  aligned  with  this  axis  so  the  decrease  in  modulus  is  linear  and  proportional  to 
the  porosity.  The  elastic  modulus  transverse  to  the  longitudinal  axis  is  more  sensitive  to  changes  in 
porosity.  The  results  for  Bert's  model  are 
for  a  hexagonal  array  of  cylindrical  voids 
and  show  a  slightly  greater  decrease  in 
modulus  with  porosity  than  the  other  two 
approaches.  The  rate  of  decrease,  however, 
is  greatest  for  the  finite  element  results. 

The  degree  of  anisotropy  predicted 
by  the  models  is  shown  by  plotting  the  ratio 
of  the  axial  modulus  to  the  transverse 
modulus.  Bert's  model  predicts  greater 
anisotropy  due  to  the  lower  transverse 
modulus  values.  The  ratio  increases  fi'om 
approximately  1.0  to  1.35  for  Zhao  et  al. 
and  the  finite  element  model,  but  reaches  a 
maximum  of  1.55  for  Bert.  These  are 
generally  lower  than  typical  ratios  reported 
for  human  compact  bone.  Cowin  [9] 
summarized  data  fi’om  four  sources  which 
yield  the  following  ratios:  1.48,  1.46,  2.66, 
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Figure  6.  Predictions  of  modulus  vs.  porosity. 
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2.16,  1.67,  1.49.  Six  values  are  included  because  two  of  the  studies  found  the  bone  to  be  orthotropic 
so  two  different  transverse  moduli  were  reported.  None  of  the  studies  included  any  microstructural 
information  about  the  specimens  tested  so  the  porosities  are  not  known.  An  important  pomt  to 
consider  as  well  is  that  the  two  analytical  approaches  assume  the  solid  phase  to  be  homogeneous  and 
isotropic.  Each  component  of  the  solid  phase  in  the  finite  element  model  is  likewise  ffeated  as 
isotropic.  The  observed  degree  of  anisotropy  is  likely  due  in  part  to  the  microstructural  architecture  of 
the  solid  phase  as  well  as  the  presence  of  voids.  The  theoretical  models  nevertheless  clearly 
demonstrate  that  a  significant  portion  of  the  anisotropy  of  Haversian  compact  bone  can  be  attributable 
to  porosity. 

Cement  line  effects.  As  mentioned  previously,  a  somewhat  mysterious  substance  known  as 
the  cement  line  separates  osteons  from  the  surrounding  interstitial  bone.  The  exact  mechanical 
properties  (modulus  and  Poisson’s  ratio)  of  the  cement  line  are  not  known  fi-om  expenmental 
testing  and  can  only  be  inferred  from  indirect  testing  and  from  the  observed  composition  of  the 
substance.  It  is  generally  believed  that  the  cement  line  is  more  compliant  [10]  than  either  the 
osteon  or  the  interstitial  bone  and  that  it  may  in  fact  be  viscous  in  nature  [11].  Moreover,  the 
precise  role  of  the  cement  line  in  determining  the  macroscopic  mechanical  behavior  of  compact 
bone  is  not  known  and  much  debated.  Accordingly,  cement  line  modulus  values  ranging  widely 
(6.0  to  0.12  GPa,  or,  50%  to  1%  of  the  osteon  modulus)  have  been  modeled  and  the  resultant 
macroscopic  properties  documented.  The  thickness  of  the  cement  line  was  also  varied  from 
1.4|j,m  to  5.0|j,m  since  it  is  reported  to  be  between  l|im  and  Sp-m  [12].  The  effects  of  these 
changes  on  macroscopic  properties  are  shown  in  Table  4  (1.4p.m  cement  line)  and  Table  5.  As 
expected,  the  cement  line  affects  the  transverse  moduli  (E22)  much  more  than  the  longitudinal 
moduli  (ki).  The  transverse  moduli  are  particularly  sensitive  to  cement  line  thickness.  Another 
interesting  observation  is  that  these  results  suggest  that  the  cement  line  modulus  is  most  likely  not 
as  low  as  0.12  GPa  (i.e.  1%  of  the  osteon  modulus),  because  the  macroscopic  transverse  modulus 
is  lower  (6.38  &  3.66  GPa)  than  generally  observed  experimentally.  Furthermore,  the  longitudinal 
modulus  results  are  essentially  the  same  as  given  by  a  rule-of-mixtures  (uniform  strain)  approach. 
The  cement  line  is  a  very  small  percentage  of  the  total  solid  phase  components  (ranging  from 
1.7%  to  8.6%),  so  the  variation  in  the  longitudinal  modulus  with  cement  line  properties  and 
thickness  is  minimal. 


Table  4  -  Effect  of  cement  line  modulus  on  macroscopic  properties 


Ec,  (GPa) 

E,i(GPa) 

E22(GPa) 

^23 

V21 

6.0  (50%)* 

12.3 

11.5 

0.31 

0.28 

3.6(30%) 

12.1 

11.3 

0.31 

0.28 

1.2(10%) 

11.9 

10.4 

0.33 

0.26 

0.12(1%) 

11.6 

6.38 

0.36 

0.16 

*  percent  of  osteon  modulus 


Table  5  -  Effect  of  cement  line  thickness  on  macroscopic  properties 


Ed 

(GPa) 

Ell  (GPa) 

E22(GPa) 

V23 

V 

21 

1.4^m 

5.0|im 

1.4|xm 

5.0(im 

1.4p,m 

5.0pm 

1.4pm 

5.0pm 

6.0  (50%)* 

12.3 

11.8 

11.5 

10.9 

0.31 

0.32 

0.28 

0.29 

3.6(30%) 

12.1 

11.6 

11.3 

10.3 

0.31 

0.33 

0.28 

0.27 

1.2(10%) 

11.9 

11.4 

10.4 

8.53 

0.33 

0.36 

0.26 

0.23 

0.12(1  %) 

11.6 

11.3 

6.38 

3.66 

0.36 

0.35 

0.16 

0.10 

*  percent  of  osteon  modulus 


Analytical  effective  medium  approaches  can  also  be  adapted  to  more  realistically  model 
the  solid  portion  of  the  bone  microstructure  as  distinct  components  (osteon,  cement  line,  «& 
interstitial  bone)  rather  than  as  a  single  phase.  The  work  of  Dasgupta  and  Bhandarkar  [13] 
considers  the  Mori-Tanaka  technique  for  the  case  of  a  fiber  surrounded  by  one  or  more 
'interphases'  embedded  in  a  matrix  material.  This  model  has  been  employed  for  the  case  of 
Haversian  cortical  bone  by  regarding  the  Haversian  canal  as  the  fiber  (modulus  approximately 
zero)  with  the  osteon  and  cement  line  as  two  interphase  materials.  Results  from  this  approach  are 
compared  with  those  of  the  finite  element  unit  cell  analysis  in  Fig.  7.  The  ratios  of  the 
longitudinal  modulus  to  the  transverse  modulus  are  plotted  vs.  cement  line  modulus  for  two 
cement  line  thicknesses.  The  two  methods  match  well  for  higher  cement  line  moduli,  but  the 
Mori-Tanaka  approach  seriously  underestimates  the  transverse  modulus  otherwise.  The  modulus 
ratios  plotted  provide  an  indication  of  the  degree  of  anisotropy  in  the  macroscopic  modulus 
predictions  and  the  effects  of  the  cement  line  upon  these.  The  longitudinal  modulus  is  only  10  to 
20  percent  higher  than  the  transverse  modulus  for  cement  line  moduli  of  50%  and  30%  of  the 
osteon  modulus.  As  the  cement  line  modulus  further  decreases,  however,  the  ratio  increases 
rather  dramatically,  particularly  for  the  thickest  (5mm)  cement  line.  Overall,  these  ratios  are 
lower  than  those  typically  reported  in  the  literature.  As  noted  previously.  Cowin's  [9]  summary 
yields  the  following  ratios;  1.48,  1.46,  2.66,  2.16,  1.67,  and  1.49. 

•At  least  two  views  are  possible  when  comparing  the  modeling  results  with  experimental 
data.  One  is  that  the  model  simply  lacks  sufficient  realism  in  its  present  form  because  it  predicts  a 
much  lower  modulus  ratio  for  most  of  the  range  of  parameters  studied.  In  this  case,  other  sources 
of  anisotropy  may  need  to  be  included  in  the  model.  Prominent  among  these  would  be  to 
consider  the  osteon  to  be  anisotropic  due  to  the  preferential  orientation  of  collagen  and  apatite  in 
the  concentric  layers  of  osteonic  lamellar  bone.  This  can  be  modeled  explicitly  within  the 
framework  of  the  present  approach  but  the  precise  degree  of  anisotropy  of  the  osteon  is  unknown. 
This  becomes  to  a  large  extent  then  another  adjustable  parameter  that  can  be  varied  over  a  wide 
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range  as  needed.  Osteon  anisotropy  has  been  investigated  theoretically  [14,15],  but  no  direct 
experimental  evidence  is  currently  available  to  firmly  establish  anisotropic  osteon  properties. 
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Figure  7.  Cement  line  effects  on  degree  of  anisotropy. 


Another  view  is  that  the  lower  cement  line  modulus  cases  in  the  present  analysis  more  accurately 
reflect  the  true  nature  of  the  cement  line  and  this  is  why  the  predictions  match  the  experimental 
data  more  closely  for  these  cases.  The  current  model  treats  the  cement  line  as  perfectly  bonded  to 
the  surrounding  bone,  which  may  therefore  overstiffen  the  model  relative  to  the  actual  physical 
situation.  The  lowest  cement  line  modulus  values  of  5%  and  1%  could  be  considered  somewhat 
unrealistic  since  they  are  20  and  100  times,  respectively,  lower  than  the  osteon  modulus,  but  the 
net  effect  taken  together  with  the  perfect  bonding  assumption  may  in  fact  combine  to  produce  a 
response  quite  close  to  the  actual  case  (because  of  slip  or  viscous  behavior).  In  either  case,  the 
overriding  need  is  for  more  definitive  experimental  data  on  the  mechanical  properties  of  the 
osteon  and  the  cement  line.  Such  information  will  not  be  easy  to  come  by,  however,  and  these 
results  demonstrate  the  utility  of  micromechanics  modeling  in  exploring  the  possible  effects  of 
microconstituent  properties  on  macroscopic  properties.  Specifically,  the  results  herein  suggest 
that  the  cement  line  is  likely  a  significant  contributor  to  Haversian  compact  bone  anisotropy.  The 
macroscopic  anisotropy  in  Haversian  bone  mechanical  properties  is  not  necessarily  due^  solely  to 
osteon  anisotropy  and  can  conceivably  be  explained  based  upon  cement  line  properties  alone. 
These  possibilities  should  therefore  remain  open  to  inquiry  until  direct  experimental  evidence 
indicates  otherwise. 
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Experimental  Testing 

The  experimental  portion  of  this  research  has  taken  several  forms  in  an  attempt  to 
characterize  macroscopic  mechanical  properties  and  their  relation  to  micro  structural  details.  The 
most  comprehensive  testing  has  been  ultrasonic  through-transmission  testing  to  measure  the 
anisotropic  elastic  moduli.  Tensile  testing  to  failure  has  also  been  conducted  to  determine  elastic 
modulus  (longitudinal  direction  only)  and  tensile  strength.  Dynamic  mechanical  testing  has  also 
been  initiated  in  order  to  measure  viscoelastic  moduli  and  damping  properties.  All  testing  has 
been  folllowed  by  extensive  microstructural  analysis  in  order  to  quantify  various  features  of  the 
bone  tissue  microstructure.  A  working  arrangement  has  been  developed  between  the  project  and 
several  laboratories  in  order  to  establish  an  experimental  protocol  for  specimen  acquistion, 
storage,  machining,  and  testing.  Fresh  equine  bone  samples  are  obtained  from  the  necropsy  lab  m 
the  Department  of  Veterinary  Physiology  and  Pharmacology  at  Texas  A&M  University.  Sections 
of  the  middle  portion  of  the  cannon  bone  were  cleaned  and  cut  to  lengths  of  approximately  2 ‘A 
inches.  The  specimens  were  kept  frozen  until  further  cutting.  They  are  then  rough  cut  into 
quarters  and  placed  in  Ringer's  solution  and  refrigerated.  A  low  speed  diamond  blade  wafering 
saw  was  used  next  to  cut  test  specimens  to  dimensions  appropriate  for  the  particular  test  method. 

Ultrasonic  testine.  The  mathematical  model  for  estimating  mechanical  moduli  by 
ultrasound  is  constructed  by  substituting  the  expression  for  a  bulk  wave  into  the  Navier  equation 
from  linear  elasticity.  The  substitution  gives  three  equations,  which  form  an  eigenvalue  problem 
with  the  terms  pv^  as  roots.  A  set  of  roots  can  be  found  after  specifying  the  direction  cosines  for 
a  given  plane  wave.  Once  the  roots  for  a  particular  set  of  direction  cosines  are  found,  the 
direction  of  displacement  is  given  by  substituting  each  of  the  roots  back  into  the  set  of  matrix 
equations.  Repeating  the  process  for  all  the  possible  pairs  of  direction  cosines  yields  the  following 
relations  between  velocities,  density(p),  and  stiffness  coefficients  (Cij): 

PVl,l“Cl1  PV2,2”C22  PV3,3“C33 

P  V2,3  “  C44  P  V3,2  “  C44 
p  Vi, 3  ~  C55  P  vl.l  “Css 
p  Vi  ,2  “Css  P  V2,1  ”  Css 


C12  =  -  Css  +  V-  C1I  -  Css  +  2p  V12.I2  V"  C22  -  Css  +  2p  Vl2,12 


Cl3  =  -  Css  +  V'  C33  -  Css  +  2p  Vl3,13  V'  Cl  1  -  Css  +  2p  Vl3,13 


C23  =  -  C44  +  -J-  C33  -  C44  +  2p  V23.23  ^/-  C22  -  C44  +  2p  V23,23 
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These  relations  are  for  the  more  general  case  of  orthotropic  material  behavior.  The  velocities  Vi,i 
are  for  longitudinal  waves  and  vg  are  for  shear  waves.  The  velocities  denoted  Vij^j  are  quasi¬ 
velocities'  measured  along  axes  45°  to  the  three  orthotropic  axes.  For  transversely  isotropic 
material  behavior,  only  one  of  these  quasi-velocities  must  be  measured.  All  of  the  velocities  in  the 
model  are  assumed  to  be  phase  velocities.  Phase  velocity  is  the  velocity  of  a  point  of  constant 
phase,  whereas  group  velocity  is  the  velocity  at  which  energy  travels  through  the  material.  The 
two  are  not  equal  for  a  dispersive  material,  but  Ashman  [16]  has  shown  compact  bone  to  be  only 
mildly  dispersive  so  group  velocity  should  be  a  reasonable  estimate  of  true  phase  velocity.  This  is 
important  for  the  present  study  since  group  velocity  is  typically  measured  by  ultrasonic 
transducers.  The  engineering  properties  can  be  obtained  by  inverting  the  stiffhess  matrix  Cy  to  get 
compliance  terms  Sy  and  performing  simple  algebra. 

The  ultrasonic  testing  protocol  is  summarized  in  Fig.  8.  The  sonic  velocities  of 
longitudinal  and  shear  waves  were  determined  for  each  specimen  using  a  through  transmission 
approach.  The  classical  technique  of  measuring  travel  time  and  specimen  thickness  was  used  to 
calculate  the  velocities.  Pulse  transducers  with  a  nominal  frequency  of  5  MHz  were  used  for  both 
longitudinal  and  shear  waves.  One  longitudinal  and  two  shear  velocities  were  measured  along 
each  of  the  three  orthogonal  axes  of  the  cube  specimen.  Following  this,  a  pair  of  oblique  faces 
were  cut  at  a  45°  angle  to  the  overall  long  axis  of  the  bone  (the  3 -axis).  Three  quasi- velocities 
were  then  measured  along  this  axis  to  complete  the  data  needed  to  calculate  the  tranversely 
isotropic  properties.  All  specimens  were  tested  wet.  A  paired  cube  immediately  adjacent  to  the 
test  cube  was  cut  and  embedded  for  microstructural  analysis. 


3  -  longitudinal 


•  5  MHz  3mm  probes 
(Panametrics  V157,  V1091) 

•  longitudinal  and  shear  waves 
generated 

•  transmission  velocities  measured 
along  all  3  main  cube  axes 

•  transversely  isotropic  material 
symmetry  assumed 

•  transmission  “pseudo-velocities” 


Figure  8.  Ultrasonic  testing  details. 
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Tensile  testing.  Specimens  for  tensile  testing  were  first  cut  into  'slabs'  measurmg 
approximately  2"x  3/8"xl/8"  on  a  low  speed  diamond  wafering  saw  under  continuous  irrigation 
with  water.  The  long  dimension  is  aligned  with  the  long  axis  of  the  bone.  These  slabs  were  then 
submitted  to  the  Testing,  Machining,  and  Repair  Facility  of  the  Texas  Engineenng  Expenment 
Station  for  machining  into  standard  'dogbone'  shaped  tensile  specimens  (see  Fig.  9).  e 
machining  is  done  at  low  speed  and  constant  irrigation  with  distilled  water.  Once  machined,  each 
specimen  is  stored  in  a  separate  vial  of  Ringer's  solution  and  refrigerated  until  tested.  Testing 
must  follow  in  less  than  90  days  or  the  specimens  should  be  refrozen.  Custom  extensometers 
were  obtained  for  measuring  the  extension  of  the  gage  section.  The  specimens  were  mounted  in 
an  Instron  servo-controlled  mechanical  testing  machine  in  articulating  grippers  to  prevent 
unintended  bending  loading.  Quasi-static  tensile  loading  was  applied  in  stroke  control  at  a  cross¬ 
head  displacement  rate  of  0.2"/minute.  A  sample  stress-strain  curve  is  presented  in  Fig.  10.  The 
elastic  modulus  is  calculated  from  the  slope  of  a  regression  line  fit  to  the  linear  portion  of  the 
curve,  and  the  tensile  strength  is  taken  to  be  the  maximum  stress  reached  during  the  test. 
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Figure  9.  Tensile  test  specimen  (all  dimensions  in  inches). 
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E2028 

Max  Stress:  134  Mpa 


Figure  10.  Typical  stress-strain  curve  from  tensile  testing. 

D^namicjnechaiiical^jMting.  Dynamic  mechanical  testing  is  commonly  employed  for 
measuring  the  viscoelastic  behavior  of  polymeric  materials.  While  there  are  many  different  tests 
available  for  characterizing  the  viscoelastic  behavior  of  materials,  such  as  creep  and  relaxation 
tests,  the  dynamic  mechanical  test  procedure  focuses  on  the  upper  boundary  of  the  frequency 
spectrum.  In  a  dynamic  mechanical  test,  an  oscillatory,  or  periodic,  strain  and  the  resulting  load 
and  stress  are  measured.  In  the  current  ivork  torsional  loading  was  applied,  in  which  case  the 
input  strain  has  the  following  form; 


y  =  rJin(cot) 

where  Yo  is  the  maximum  strain  amplitude  and  O)  is  the  angular  frequency  in  radians  per  second. 
For  linear  viscoelastic  behavior,  when  equilibrium  is  reached  the  shear  strain  and  shear  stress  both 
vary  sinusoidally,  but  the  stress  lags  beWnd  the  strain.  The  resulting  output  for  strain-controlled 
loading  is  a  shear  stress  of  the  following  form; 


r  =  ToSin(ctyt  +  S) 

where  to  is  the  amplitude  of  the  shear  stress.  The  phase  lag  between  the  shear  stress  and  applied 
shear  strain  is  given  by  5.  The  phase  lag  ranges  from  90°  for  a  purely  viscous  material,  to  0°  for  a 
purely  elastic  material.  In  a  dynamic  mechanical  test,  t  and  y  are  related  by. 
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t=y  JG*  =  /„[G'sin(fyO  +  G"cos(fyO] 

where  G*  is  the  complex  modulus,  and  G'  and  G"  are  defined  as  the  storage  modulus  and  loss 
modulus,  respectively,  and  are  given  by; 


G  =  — cos<5 

Yo 


and 


G"=^s\n5 

Yo 


The  loss  factor,  or  tan  S,  can  then  be  expressed  as  a  function  of  G  and  G": 


o  G" 
tano  =  — 
G 


It  is  then  obvious  that  as  the  phase  angle  approaches  zero,  the  storage  modulus  approaches  the 
elastic  modulus  and  the  loss  modulus  goes  to  zero.  In  dynamic  mechanical  analysis  with  torsional 
loading,  a  specified  sinusoidal  rotation  is  applied  to  one  end  of  the  specimen.  The  test  machine 
measures  the  torque  at  the  other  end  of  the  specimen,  as  shown  in  Fig.  11,  and  the  phase 
difference  between  the  stress  and  strain. 

Dynamic  mechanical  testing  was  conducted  on  a  Rheometrics  RDS  II  located  in  the 
Polymer  Technology  Center  at  Texas  A&M  University.  Specimens  for  dynamic  mechanical 
testing  consisted  of  slabs  approximately  2"x  3/8"xl/8"  in  size.  Preliminary  tests  were  conducted 
to  determine  an  appropriate  strain  amplitude  and  frequency  range  for  linear  viscoelastic  behavior. 
Specimens  were  stored  in  Ringer's  solution  and  refrigerated  until  just  prior  to  testing.  They  were 
then  allowed  to  reach  room  temperature  but  maintained  moist.  Immediately  before  mounting  in 
the  testing  machine,  each  specimen  was  coated  with  a  thin  layer  of  high  vacuum  silicone  grease  to 
prevent  moisture  loss  during  testing.  The  torsional  strain  amplitude  was  0.015%  and  the 
frequency  was  varied  over  three  decades  from  0.1  to  100  rad/s.  Each  test  took  roughly  10 
minutes.  Results  are  reported  as  storage  modulus  (G'),  loss  modulus  (G"),  and  loss  factor  (tan  6) 
as  a  function  of  frequency.  A  typical  set  of  plots  is  shown  in  Fig.  12. 
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Microstructural  analysis.  Specimens  from  all  three  testing  types  were  analyzed  for 
microstructural  details.  Quantitative  characterization  of  the  microstructure  can  be  grouped  into 
two  general  categories;  (1)  physical  properties,  such  as  densities  and  mineralization;  and  (2) 
image  analysis  of  microradiographs  of  transverse  sections. 
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The  first  quantity  measured  was  the  bulk  density  -  the  average  of  the  bone  composition, 
water  content,  and  void  spaces.  A  method  based  upon  Archimedes  Principle  allows  the  wet 
density  to  be  determined  by  the  following  relation: 


Phone 


moist_ 


[IV  -W  J 

moist  ^'susp'' 


where,  pbone  is  the  bulk  wet  density  of  the  specimen,  Wmoist  is  the  weight  of  the  moist  specimen, 
Wsusp  is  the  weight  of  the  specimen  when  it  is  suspended  with  wire  in  distilled  water,  and  Pwater  is 
the  density  of  distilled  water.  An  analytic  balance  was  used  to  weigh  the  specimens  with  a 
precision  of  five  decimal  places  when  measuring  grams.  The  dry  density  was  determined  next. 
The  volume  of  the  section  was  calculalted  the  wet  density  calculations,  since  the  density  and  wet 
weight  are  known,  and  also  measured  mechanically  using  a  digital  caliper.  The  specimens  were 
placed  in  pyrex  beakers  and  dried  in  a  convection  oven  at  60°Celsius  for  about  two  days.  The 
specimens  were  then  weighed  again  using  the  analytic  balance  by  weighing  the  specimen  and 
beaker,  together,  then  subtracting  the  weight  of  the  beaker,  thus  giving  the  weight  of  the  dried 
specimen.  The  dry  density  was  then  be  computed  as  the  ratio  of  the  dry  weight  to  the  volume  of 
the  specimen: 

_  ^dry 
Pdry  ~  y 


where,  Wdry  is  the  dry  weight  of  the  specimen,  while  V  is  the  volume  of  the  specimen. 


Another  microstructural  parameter  measured  was  the  mineral  content  for  each  specimen. 
In  order  to  quantify  the  mineral  content,  the  specimen  must  first  be  ashed.  This  process  involved 
cooking  the  specimens  (in  their  pyrex  beaker  with  crucible  lids  covering)  at  600°C  for  24  hours  in 
order  to  bum  off  all  other  material,  thus  reducing  the  bone  to  its  mineral  phase.  This  remaining 
“ash”  is  composed  of  a  combination  of  calcium  crystals,  phosphorus,  and  other  trace  materials. 
The  specimens  were  then  weighed  again  on  the  analytic  balance  after  the  oven  was  allowed  to 
cool  to  room  temperature.  Cooling  the  specimens  is  cmcial  because  the  mass  will  fluctuate  until 
the  specimens  reach  room  temperature.  The  %  mineralization  is  then  calculated  as  the  ratio  of 
the  weight  of  the  remaining  ash  material  to  the  dry  weight  of  the  whole  specimen: 


VoMineralization  = 


''ash 


*100 


where,'Wdry  is  the  dry  weight  of  the  specimen,  and  Wash  is  the  weight  of  the  ash  material  of  the 
specimen.  This  quantity  is  also  called  the  “ash  weight”  in  the  literature,  or  “ash  weight  %”,  even 
though  the  units  are  a  simple  percentage. 


Microradiography  and  image  analysis  were  conducted  through  subcontract  arrangements 
with  the  University  of  Texas  Medical  Branch  (UTMB)  at  Galveston.  Samples  from  each  test 
specimen  were  cut  and  imbedded  in  Buehler  Castolite  Resin.  Next,  a  100  micron  thick  slice  was 
cut  using  an  Isomet  saw.  A  microradiograph  was  then  taken  using  a  10  kilovolt,  20  millamp 
Sorensen  x-ray  system,  with  an  exposure  time  of  approximately  10  minutes.  Fig.  12  shows  a 
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sample  microradiagraph  taken  of  a  section.  The  porous  portions  of  the  bone  specimen 
microradiographed  are  indicated  by  the  darkest  regions.  Conversely,  the  most  highly  mineralized 
portions  of  the  bone  are  illustrated  by  the  regions  with  the  lightest  colors.  The  porosity  of  the 
bone  tissue  in  each  microradiograph  image  was  using  Optimas  4.2  image  analysis  software. 
Threshold  levels  were  adjusted  manually  until  all  'dark'  regions  were  highlighted  and  then  the 
percent  area  of  the  highlighted  pixels  was  calculated.  Each  field  is  1  mm^  in  area.  The  number  of 
fields  analyzed  for  each  specimen  varied  with  specimen  size.  The  number  of  fields  analyzed  was  9 
for  the  specimens  from  ultrasonic  testing,  5  from  the  tensile  specimens,  and  8  from  the  specimens 
from  the  dynamic  mechanical  testing.  Results  were  then  averaged  for  all  fields  analyzed.  Other 
microstructural  parameters  of  interest  were  the  extent  of  cement  line  material  present  and  the  area 
fraction  of  the  cross-section  occupied  by  Haversian  systems  (osteons).  These  were  not  calculated 
for  all  cases  because  the  procedures  are  much  more  labor  intensive.  In  this  case,  the  outer 
boundaries  of  each  osteon  in  the  field  was  traced  manually  using  the  'mouse'  on  the  PC.  The  total 
length  of  the  lines  traced  was  then  calculated  by  the  image  analysis  software.  The  length  was  then 
normalized  by  dividing  by  the  total  area  of  the  field  (1  mm')  to  give  a  cement  line  'density'  (in 
units  of  mm/mm').  The  area  enclosed  by  each  closed  curve  was  also  calculated  and  the  percent 
osteon  area  determined  as  this  quantity  divided  by  the  total  area  of  the  field  (1  mm  ). 


Figure  12.  Sample  microradiograph  of  bone  section  (63X). 


Properties  vs.  Microstructure 


In  order  to  gain  insight  into  trends  and  effects,  results  from  all  three  modes  of  mechanical 
testing  were  analyzed  and  plotted  versus  various  microstructural  parameters. 

Anisotrnnic  moduli.  The  main  results  from  the  ultrasonic  tests  were  a  M  set  of 
anisotropic  moduU  (and  Poisson's  ratios)  for  each  specimen,  assuming  transversely  isotropic 
material  symmetry.  The  variation  in  longitudinal  and  transverse  elastic  moduli  with  porosity  is 
depicted  in  Fig  14  The  trends  are  as  expected,  with  moduli  generally  decreasing  with  increasing 
porosity.  No  data  has  been  previously  reported  for  variation  in  transverse  modulus  with 
microstructure,  so  the  longitudinal  modulus  alone  is  compared  with  other  experimental  results  m 
Fig.  15.  Note  that  the  data  of  Schaffler  and  Burr  [7]  has  been  divided  into  osteonal  and  ple^afom 
tissue  types.  The  dashed  line  representing  a  linear  regression  is  fit  only  to  the  osteonal  results  for 
more  direct  comparison  with  the  present  results.  Data  from  a  recent  study  by  McCalden,  etal. 
[17]  is  also  included  (but  only  plotted  for  the  same  range  of  porosities  as  the  current  work). 
Further  these  results  were  generated  by  manuaUy  cross-plotting  modulus  vs.  porosity  from  ^o 
separate  graphs  of  modulus  vs.  age  and  porosity  vs.  age  as  presented  m  the  ongmal  paper.  Ihe 
regression  line  for  the  present  study  has  a  slope  between  the  other  two  data  sets.  An  unexpected 
feature  of  the  present  results  (Fig.  14)  is  that  the  transverse  modulus  appears  to  be^ 
to  variations  in  porosity  than  the  longitudinal  modulus.  This  is  contradictory  to  w  at  w^ou  e 
expected  based  upon  an  idealized  theoretical  treatment  of  porosity  as  parallel  hollow  cylmdnca 
voids.  The  theoretical  results  presented  previously  are  compared  to  the  current  expenmental 
results  in  Fig.  16.  The  theoretical  results  show  that  both  the  longitudinal  (E33)  and  transverse 
(E22)  moduli  do  indeed  decrease  with  porosity,  but  the  ratio  E33/E22  increases  with  PO^osity 
highlighting  the  greater  sensitivity  in  E22.  This  confounding  finding  suggests  that  such  a  simplified 
view  of  the  architecture  of  Haversian  canals  (and  other  void  spaces  as  well)  may  be  inadequate. 

The  variation  in  shear  moduli  with  porosity  is  summarized  in  Fig.  17.  Again,  both  sets  of 
moduli  generally  decrease  with  increasing  porosity,  as  would  be  expected.  These  results  also 
show  however,  that  the  sensitivity  to  variations  in  porosity  is  not  the  same  for  the  out-ot-p  ane 
shear  moduli  (G13 ,  G23 )  and  the  'in-plane'  shear  moduli  (Gn  )•  Specifically,  the  larger  (negative) 
slope  for  the  out-of-plane  moduli  indicates  greater  sensitivity  to  porosity.  Plots  showng  t  e 
variation  in  elastic  moduli  and  shear  moduli  with  percent  ash  weight  are  given  in  Figs.  18  and  y, 
respectively.  All  moduli  generally  increase  as  ash  weight  increases,  but  the  longitudinal  modulus 
is  more  sensitive  than  the  the  transverse  modulus  and  the  out-of-plane  shear  modulus  is  more 
sensitive  than  the  in-plane  shear  modulus.  Similar  plots  are  included  for  the  vanation  in  rno  u  1 
with  osteon  area  fraction  and  cement  line  perimeter  in  Figs.  20-23.  These  results  are  consistent 
with  previous  qualitative  observations  that  secondary  osteons  tend  to  decrease  the  modulus  an 
strength  of  compact  bone  tissue  [18-20].  The  continual  remodeling  that  produces  secondary 
osteons  would  be  expected  to  increase  both  the  osteon  area  fraction  as  we  as  t  e  cemen 
perimeter.  Once  again,  the  longitudinal  modulus  and  out-of-plane  shear  modulus  show  ^eater 
sensitivity  to  microstructural  variation,  increasing  osteon  area  fraction  and  cement  line  perimeter 

in  this  case. 
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Longitudinal  Modulus  -  E33  (GPa)  Elastic  Moduli  (GPa) 


Porosity  (%) 

Figure  14.  Elastic  moduli  vs.  porosity. 
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Figure  15.  Comparing  longitudinal  modulus  variation  with  others, 
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Figure  16.  Experiment  vs.  theory  for  normalized  moduli. 
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Figure  17.  Shear  moduli  vs.  porosity. 
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Figure  22.  Elastic  moduli  vs.  cement  line  perimeter. 
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Figure  23.  Shear  moduli  vs.  cement  line  perimeter. 
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Tensile  testin2.  A  total  of  12  'dogbone'  tensile  specimens  have  been  tested  thus  far  in  the 
project.  An  advantage  to  tensile  testing  is  the  ability  to  measure  strength  and  related  pwameters 
in  addition  to  elastic  modulus.  The  longitudinal  elastic  moduli  and  tensile  strengths  ^e  plotted  vs^ 
porosity  in  Figs.  24  and  25,  respectively.  Both  modulus  and  strength  decrease  with  porosity  but 
the  negative  slope  for  the  modulus  seems  to  be  steeper  than  reported  above  using  ultrasomc 
methods  for  measuring  moduli.  Furthermore,  the  correlation  for  moduus  is  rat  er  we  as 
indicated  by  the  dashed  line,  but  is  much  better  when  the  results  for  the  two  lowest  porosity 
specimens  are  not  included.  Additional  microstructural  analysis  is  also  being  conducted  for  the 
tensile  test  specimens.  Twelve  more  specimens  have  already  been  machined  and  are  awaiting 
testing,  so  the  appropriateness  of  these  two  data  points  should  become  more  apparent  as  more 
results  are  obtained.  In  order  to  estimate  the  degree  of  uniformity,  or  lack  thereof,  in  the 
mineralization  levels  within  the  tissue  microstructure,  the  mineralization  distnbution  is  being 
determined  from  the  microradiograph  images.  The  Optimas  4.2  image  analysis  software  is  again 
employed  in  this  effort.  During  the  microradiographing  procedure,  an  alummum  step  wedge  is 
also  exposed  with  each  slice  of  bone.  This  provides  gray  scale  levels  for  de  mng  i  erent 
degrees  of  mineralization  as  'lowly',  'moderately',  or  'highly'  mineralized  bone  tissue^  This  data, 
along  with  cement  line  perimeter  results,  should  provide  considerable  ^sight  into  the 
microstructural  features  and  mechanisms  responsible  for  strengthening  and  stiffemng  Haversian 
compact  bone  tissue.  Because  these  procedures  are  fairly  labor  intensive,  they  are  currently  still  in 
progress  and  should  be  completed  in  the  near  future. 


Porosity  (%) 


Figure  24.  Modulus  variation  with  porosity. 
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Dynamic  mechanical  properties.  The  representative  frequency  sweep  presented 
previsously  in  Fig.  12  indicates  that  the  storage  (G')  and  loss  (G")  moduli  only  change  slightly 
over  the  range  of  frequencies  tested.  The  loss  factor  (tan  5)  decreases  over  the  first  two  decades 
and  then  levels  off  or  increases  slightly.  Results  from  the  frequency  sweep  for  each  specimen 
were  thus  summarized  in  the  form  of  average  values  for  the  purposes  of  correlating  with 
microstructural  parameters.  The  storage  and  loss  moduli  are  plotted  vs.  wet  and  diy  density  in 
Figs.  26  and  27,  respectively.  The  storage  modulus  increases  with  density,  as  would  be  expected 
but  the  loss  modulus  somewhat  surprisingly  remains  essentially  constant  over  the  range  of 
densities  tested.  Fig.  28  shows  a  similar  trend  for  moduli  vs.  %  mineralization.  Since  the  purpose 
of  these  tests  is  to  seek  to  detect  the  influence  of  the  cement  line  on  dynamic  properties,  the  loss 
factor  (which  is  the  ratio  of  G"  to  G’)  is  plotted  vs.  cement  line  perimeter  in  Fig.  29.  The  results 
reveal  a  fairly  consistent,  although  mild,  decrease  in  loss  factor  with  increasing  cement  Imes^  Itus 
is  also  opposite  the  expected  trend,  at  least  under  the  hypothesis  that  the  cement  line  would  tend 
to  increase  damping  within  the  material.  The  storage  and  loss  moduli  are  plotted  vs.  /o  secondary 
Haversian  area  (SHA)  in  Figs.  30  and  31.  Fig.  30  shows  that  many  samples  have  very  low 
%SHA,  which  suggests  that  they  were  mainly  primary  bone  and  lacked  significant  Haversian 
microstructure.  Fig.  3 1  shows  the  same  results  but  for  %SHA  greater  than  25%.  In  this  case  the 
storage  modulus  does  in  fact  decrease  with  increasing  %SHA  as  has  been  observed  by  others  (and 
was  also  confirmed  in  our  ultrasonic  tests  reported  previously  herein).  Thus,  more  conclusive 
examination  of  the  role  of  the  cement  line  in  damping  awaits  testing  of  more  samples  with 
predominately  Haversian  microstructure. 
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Figure  26.  Variation  in  viscoelastic  moduli  with  wet  density. 


Figure  27.  Variation  in  viscoelastic  moduli  with  dry  density. 
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Figure  28.  Variation  in  viscoelastic  moduli  with  ash  weight  (%  mineralization). 


Figure  29.  Variation  in  loss  factor  with  cement  line  perimeter. 
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Figure  30.  Variation  in  loss  factor  with  %  secondary  Haversian  area  (SHA). 


Figure 


3 1 .  Variation  in  loss  factor  with  %  SHA  for  predominantly  Haversian  bone  only. 


Loss  Modulus  (GPa) 


Dissemination  of  Results 
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publications  as  listed  below.  Three  full-length  refereed  journal  articles  are  also  currently  in 
preparation. 

Ayers,  A.  K.,  and  Hogan,  H.  A.,  "Dynamic  Mechanical  Properties  of  Equine  Compact  Bone  and 
Their  Dependence  Upon  Tissue  Density,"  accepted,  1995  Bioengineering  Conference,  June 

28-July  2,  1995,  Beaver  Creek,  CO.  .. 

Hogan,  H.  A.,  Guna  K.  S .  Ahern,  J,  C.,  and  Simmons,  D,  J,,  "Charact^g  Oaeon^ 
Microstructure  and  Its  Influence  Upon  the  Anisotropic  Moduli  of  Equme  Cortical  Bone, 
accented.  1995  Bioengineering  Conference,  June  28-July  2,  1995,  Beaver  Creek,  CO. 

Ayers  A  K  and  Hogan,  H.  A.,  "Correlation  of  Viscoelastic  Properties  and  Microstructure  for 
Equine  Cortical  Bone  Tissue,"  Abstracts  of  the  Thirteenth  Annual  Houston  Conference  on 

Biomedical  Engineering  Research,  February  16-17,  1995,  Houston,  TX,p.  97.  ,  .  .  , 

Hogan,  H.  A.,  Ahem,  J.  C.,  and  Simmons,  D.  J.,  "Variations  in  the  Anisotropic  Mech^cal 
Properties  of  Compact  Bone  with  Tissue  Micro  structure,"  Advances  in  Bioengineenng  1994, 
BED  Vol.  28,  1994  ASME  Winter  Annual  Meeting,  Nov.  6-11,  1994,  Chicago,  IL,  pp.  229- 

DeFrlt,  R.  J.,  Ahem,  J.  C.,  and  Hogan,  H.  A.,  "The  Anisotropic  Mechanical  Properties  of 
Compact  Bone  and  Their  Correlation  with  Tissue  Microstmcture,"  Recent  Ad^nces..._i|i 
Engineering  Science,  proceedings  of  the  Society  of  Engineering  Science  31*  Annual  Techmcal 

Meeting,  October  10-12,  1994,  College  Station,  TX,  pp.  255.  ,  •  r?i  *• 

Ahem,  J.  C.,  Hogan,  H.  A.,  and  Simmons,  D.  J.,  "Relations  Between  the  Ultrasomc  Elastic 
Moduli  of  Compact  Bone  and  Tissue  Microstmcture,"  Biomedical  Engineering  -  Recent 
Developments,  proceedings  of  the  Thirteenth  Southern  Biomedical  Engineenng  Conference, 

April  16-17,  1994,  Washington,  DC,  pp.  681-684. 

DeFrese,  R.  J.,  and  Hogan,  H.  A.,  "The  Effect  of  Imperfect  Cement  Lme-Osteon  Bonding  on  the 
Properties  of  Cortical  Bone  Tissue,"  Biomedical  Engineering  -  Recent  Developi^tits, 
proceedings  of  the  Thirteenth  Southern  Biomedical  Engineering  Conference,  Apnl  16-17, 

1994,  Washington,  DC,  pp.  685-688.  .  ,  ••  fTj 

Ahem,  J.  C.,  and  Hogan,  H.  A.,  "A  Study  of  Dispersion  in  the  Ultrasonic  Velocities  of  Havesiari 
Compact’  Bone,"  Abstracts  of  the  Twelfth  Annual  Houston  Conference  on  Biomedical 

Engineering  Research,  Febmaiy  10-11,  1994,  Houston,  TX^  p.  93. 

DeFrese,  R.  J.,  and  Hogan,  H.  A.,  "Micromechanics  Modeling  of  Compact  Bone:  Interphase 
Effects,"  Abstracts  of  the  Twelfth  Annual  Houston  Conference  on  Biomedical  Engineering 

Research, Febraary  10-11,  1994,  Houston,  TX,  p.  94.  ,  .  j  i-  e 

Hogan,  H.  A.,  De  Frese,  R.  J.,  and  Simmons,  D.  J.,  "Composite  Micromechamcs  Modelmg  of 
Porosity  and  Interphase  Effects  in  Compact  Bone  Tissue,"  presentation  to  Seminar  m 

Mechanics  and  Materials  -  MEEN  681  c\a.ss,  January  27,  1994.  ...  ■  c 

Hogan,  H.  A.,  "Some  Biomimetics  Problems  in  Biomechanics,"  Invited  Seminar,  joint  meeting  ot 
A&M  Student  Chapter  of  the  Biomedical  Engineering  Society  and  Bioengineering  Graduate 
Seminar  Class,  Texas  A&M  University,  September  16,  1993,  College  Station. 


30 


Hogan  H  A  and  DeFrese,  R.  J.,  "Theoretical  Modeling  of  the  Effects  of  Porosity  on  the 
Mechanical  Properties  of  Compact  Bone  Tissue,"  TSB  ’93  Resumes/Abstracts,  Vol.  I, 
International  Society  of  Biomechanics  XIV*  Congress,  July  4-8,  1993,  Paris,  FRANCE,  pp. 


582-583 

DeFrese  R  J.  and  Hogan,  H.  A.,  "A  Micromechanics  Study  of  the  Influence  of  the  Cement  Line 
on  the  Anisotropy  of  Compact  Bone  Tissue,"  1993  Bioengineering  Conference,  BED  Vol.  24, 
N.  A  Langrana,  M.  H.  Freidman,  and  E.  S.  Grood,  eds.,  1993  ASME/AICHE/ASCE  Summer 
Bioengineering  Conference,  June  25-29,  1993,  Breckenridge,  CO,  pp.  238-241. 

Hogan,  H.  A.  and  DeFrese,  R.  J.,  "A  Comparison  of  Methods  for  Including  Porosity  m 
Micromechanics  Modeling  of  Compact  Bone  Tissue,"  Advances  in  Experimental  Mechanics 
andBiomimetics.  AD  Vol.  29,  AMD  Vol.  146,  1992  ASME  Winter  Annual  Meeting,  Nov.  8- 

13, 1992,  Anaheim,  California,  pp.  111-121.  .  ,  -r 

DeFrese,  R.  J.  and  Hogan,  H.  A.,  "The  Influence  of  Packing  Geometry  and  Fiber/Matnx  Interface 
Conditions  on  Composite  Micromechanics  Modeling  of  Compact  Bone  Tissue,  1992  ASME 
Applied  Mechanics,  Materials  and  Aerospace  Summer  Meeting,  Apnl  28  -May  1,  1992, 
Scottsdale,  Arizona. 


Graduate  stndpnt  research.  The  research  conducted  as  part  of  this  project  has 
constituted  graduate  thesis  work  for  the  following  three  graduate  students.  Please  note  that  the 
work  of  John  C.  Ahem  was  also  partially  supported  by  a  grant  from  the  NIH  (1  R03  RR06856- 
01A2). 

Raymond  J.  DeFrese,  "Modeling  and  Characterization  of  Haversian  Cortical  Bone  as  a  Fiber 
Reinforced  Composite  Material,"  Ph.D.,  to  graduate  December  1995. 

Andrew  K.  Ayers,  "Correlations  of  Mechanical  Viscoelastic  Properties  to  Microstmcture  of 

Equine  Cortical,"  M.S.,  August  1995.  r/-'  •  m 

John  C.  Ahem,  "Correlations  Between  the  Anisotropic  Mechanical  Properties  of  Cortical  Bone 

Tissue  and  Microstmcture,"  M.S.,  August  1994. 


31 


SUMMARY 


The  strategy  employed  in  the  research  conducted  under  this  grant  was  to  first  develop  a 
valid  theoretical  tool  for  studying  the  micromechanics  of  compact  bone  tissue,  and  then  use  this 
tool  to  study  structure/property  relations  for  potential  biomimicking  applications.  Thus,  the  early 
phases  of  the  work  focused  upon  developing  and  evaluating  fimte  element  umt  cell 
micromechanics  approaches.  Mori-Tanaka  effective  medium  methods  were  also  examine  .  ue 
to  lack  of  available  results  in  the  literature  fi-om  other  sources,  however,  focus  soon  shifted  to 
generating  experimental  data  for  verifying  and  validating  the  theoretical  models.  The  correlations 
established  from  ultrasonic  testing  suggest  that  simple  2-D  approximations  are  likely  inadequate 
for  micromechanics  modeling  of  Haversian  compact  bone  tissue.  In  particularly,  the  dependence 
of  anisotropic  moduli  upon  porosity  is  opposite  the  trend  expected  from  a  parallel  holbw  cylinder 
arrangement  of  voids  (representing  the  Haversian  canal).  Thus,  future  theoretical  efforts  should 
be  directed  at  including  three-dimensional  effects,  either  directly  through  3-D  modeling  or 
indirectly  through  modified  or  'quasi'  2-D  approaches.  These  findings  are  specifically  related  to 
predictions  of  ela.stic  moduli  and  not  necessarily  true  for  strength  or  dynamic  properties. 
Correlations  between  tensile  strength  and  microstructure  have  been  initiated  during  the  Project 
period,  and  are  continuing  currently  in  order  to  include  a  large  enough  sample  size  for  meaningful 
conclusions.  More  effort  than  originally  expected  has  been  directed  dunng  this  project  toward 
examining  dynamic  viscoelastic  behavior.  Because  the  cement  line  present  in  Haversian  compact 
bone  is  similar  to  the  compliant  interface  of  brittle  composites,  dynamic  testing  has  been  imtiated 
under  the  assumption  that  the  cement  line  will  have  greater  influence  upon  dynamic  behawor  than 
upon  static  elastic  moduli.  The  results  reported  herein  have  established  a  viable  method  for  such 
tests.  The  correlations  between  loss  factor  and  microstructure  are  ambiguous  because  of  a  large 
number  of  specimens  that  did  not  possess  exclusively  Haversian  type  tissue  microstructure. 
Further  work  is  needed  in  this  area,  both  experimentally  and  theoretically.  Expenmentally,  torsion 
testing  will  be  combined  and  compared  with  bending  testing  as  the  mode  of  loading,  and  a  l^ger 
number  of  Haversian  bone  specimens  will  be  obtained.  Both  equine  and  human  (and  perhaps 
bovine)  bone  will  be  compared  as  well.  In  summary,  much  process  has  been  made  over  the 
project  period  in  defining  the  important  features  needed  in  a  valid  micromecharacs  model  tor 
compact  bone  tissue  and  establishing  correlations  between  macroscopic  behavior  and  quantitative 
microstructural  parameters.  With  further  effort  focusing  upon  dynamic  behavior  and  matena^ 
strength,  more  insight  should  be  gained  into  the  unique  role  of  the  cement  line  interphase  matenal 
and  its  significance  for  potential  biomimicking  applications. 
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